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Abstract 

We study the discrete-time quantum walk-based search for a marked vertex on a graph. By 
considering various structures in which not all vertices are equivalent, we investigate the relationship 



between the successful search probability and the position of the marked vertex, in particular its 
' centrality. We find that the maximum value of the search probability does not necessarily increase 

in 

as the marked vertex becomes more central and we investigate an interesting relationship between 

r- 1 ! 

the frequency of the successful search probability and the centrality of the marked vertex. 
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I. INTRODUCTION 



Quantum walks are the quantum analogue of classical random walks. Rather than step- 
ping with a certain probability between adjacent vertices of a graph, a quantum walker is 
characterised by a set of probability amplitudes associated with vertices of the graph [1( • The 
strikingly different behaviour of quantum walks from their classical counterparts has already 
been harnessed in the formulation of quantum walk-based algorithms which can outperform 
corresponding classical algorithms [2- 4] . While the continued interest in quantum walks is 
largely due to these algorithmic applications in the context of quantum computation, the 
quantum walk also forms a powerful and flexible model of the evolution of a coherent (or 

fin 

partially decoherent) quantum system |5j-|9|]. Since analytical techniques are currently being 
developed to analyse quantum walks, their application to diverse problems in physics will 
likely become increasingly common. Given also the usefulness of classical random walks in 
studying transport on complex structures, it is interesting from a physical standpoint to 
continue to characterise quantum walks on graphs. 

Searching is one of the major problems in computer science and a large amount of research 
in the field of theoretical quantum computation has been in the development of general 
algorithms for fast searching of databases. Quantum search algorithms were first introduced 
by Grover to search an unsorted database [10, LL3J and later extended to quantum walk- 



based search algorithms for specific database topologies in both the discrete-time |12Ml4| 
and continuous-time 15[ cases. These studies focused on highly symmetric structures such 
as hypercubic lattices, complete and bipartite graphs, and found that the topology of the 
database was crucial in determining the efficiency of the search. An important difference 
between the discrete and continuous-time formulations of quantum walks is the extra "coin" 

n 

degrees of freedom in the discrete-time case. Ambainis et al. [13[ demonstrated that discrete- 
time search can achieve full quadratic speedup relative to classical search for hypercubic 
lattices in > 3 dimensions and it outperforms continuous-time search for lattices in two 



spatial dimensions. As shown by Childs &; Goldstone 



15], continuous-time search on the 



hypercube only achieves the full quadratic speedup for d- dimensional lattices where d > 4. 



A recent paper by Agliari et al. 16] considers continuous-time quantum walk-based search 
on fractals and thus represents the first effort to characterise the search procedure on struc- 
tures that are not vertex-transitive. An interesting phenomenon which arises when studying 



structures where not all vertices are equivalent is that the successful search probability de- 
pends on the location of the marked vertex. In their study of quantum search on Cayley 
trees, T-fractals and dual Sierpihski gaskets, Agliari et al. assumed that a peripheral vertex 
would be more difficult to find than a more central vertex, i.e. the maximum success prob- 
ability for a central vertex would be greater than for a peripheral vertex. In this work, we 
analyse this idea in more detail by studying how the maximum success probability varies 
with the centrality of the marked vertex. We find that in some simple cases, the maximum 
success probability does indeed increase with increasing centrality. However we show that, 
in general, such a relationship does not hold. 

The efficiency of quantum walk-based search relative to classical search is not only deter- 
mined by the maximum success probability, but also the time taken to reach the maximum. 
We therefore analyse the lowest frequency of the success probability as an indicator of the 
time complexity of the search. Our results suggest that this frequency is correlated with the 
centrality for a larger class of graphs than the maximum success probability and we discuss 
exceptions in terms of local structure of these graphs. 

In this article, we study discrete-time quantum walk-based search on non- vertex-transitive 
structures, and show that the frequencies present in the success probability are determined 
by the global structure of the graph as well as the centrality and local structure surrounding 
the marked vertex. To the best of our knowledge, our derivation in the Appendix contains 
the first analytical solution for discrete-time quantum walk-based search on a finite line 
with two reflecting boundaries, which provides the characteristic frequencies of the quantum 
walk-based search probability. 

The article is organised as follows. Sec. HT1 provides an introduction to quantum walks and 
quantum walk-based search. Sec. [Till describes the measure of centrality used. In Sec. [TV] 
we describe the structures considered and in Sec. [V] |VI| IVHI we present our analytical and 
numerical results. Finally, Sec. IVI III contains discussions and conclusions. In the Appendix 
we provide details of our analytical calculations. 

II. QUANTUM WALK-BASED SEARCH ON GRAPHS 

Let G(V,E) be an undirected graph with vertex set V = {v 1; v 2 , v 3 , . . .} and edge set 
E = {(vi, Vj), (vk, Vi), . . .} consisting of unordered pairs of connected vertices. If there are d 
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edges incident on a vertex Vi, we say that V{ has degree d. As described in [17], |18|, T-L-p is 
defined as the position Hilbert space, which is spanned by an orthonormal basis of vertex 
states {\vi) : G V}. For a graph of maximum degree d, He is defined as the <i-dimensional 
coin Hilbert space spanned by the orthonormal basis of coin states {|cj) : i = l,...,d}, 
representing the outgoing edges at a vertex t>j. The discrete-time quantum walk considered 
here takes place on the subnodes of the graph, which are represented by product states of 
the form \v) £g> |c) = \v,c) G H-p <g> He- Note that if G is not d-regular, then there are 
vertices of degree di < d. In this case the states {\vi, c) : c > di} do not physically represent 
subnodes of the graph G and are not occupied at any stage of a quantum walk on G. 

Let one step of the discrete-time quantum walk on the graph be the application of the 
unitary time-evolution operator U — S ■ (1 <g> C), where S is the shift operator and C is the 
coin operator. S acts on the extended position space %p <S> He as, 

S\vi,Cj) = \vj,Ci), (1) 

where \vi,Cj) is the subnode state corresponding to the edge (v^Vj) at the vertex Vi. The 
coin operator C at a vertex V{ of degree dj can be represented by a di x di matrix, which 
mixes the probability amplitudes of the subnode states of t>j. We mainly consider symmetric 
coin matrices, so that the ordering of the subnodes at a particular vertex is unimportant. 
In the Appendix, when we reduce the quantum walk on a Cayley tree to a one-dimensional 
walk, we need a biased coin. The labelling will then become important and will be made 
explicit. 



12| for the discrete-time quantum 



We follow the procedure introduced by Shenvi et al. 
walk-based search for a marked item. The quantum walker initially has equal probability 
to be found at each vertex, i.e. the state |\&o) is an equal superposition of all vertex states 
\vi) G Hp. The probability amplitude at each vertex is then divided equally between all 
subnodes. It should be noted that for a graph which is not degree-regular this is not the 
same as an equal superposition of all subnode states. Formally, the initial state is given by 

1 n <k 1 

w^EE^Ki). (2) 



Now consider a subset M C V of marked vertices. The marking is intended to represent a 
"quantum oracle" and is implemented as a perturbation to the coin operator at the marked 
vertices. A precise description of an oracle in this context can be found in 12j. The 



quantum search procedure proceeds via the repeated application of the perturbed time- 
evolution operator, U' = S ■ (1 <8> C) where the coin operator at vertex Vi is given by 

_ , Smn + 2/di ,Vi^M 
-5 mn , Vi e M 

for m,n = 1, . . . ,di. 

The coin operator above for ^ M will be referred to as the Grover coin. For all examples 
in this paper, the set M will contain only 1 vertex. The success probability, P s (t), is defined 
as the probability of finding the quantum walker at the marked vertex v m at time t. This is 
given by 

P,(t) := \(v m \*(t))\ 2 = \(v m \(U'Y\^ )\ 2 where v m G M. (4) 

The time averaged success probability is denoted by (P s ). As shown in Eq. EJ P s (t) is de- 
termined by unitary time-evolution from the initial state. Reversibility of unitary processes 
implies that P s {t) does not converge for large t but instead oscillates, which allows us to 
define the search frequency u s as the lowest frequency present in the success probability. On 
graphs for which quantum walks are not exactly periodic, this is computed from the power 
spectrum of P s (t) by selecting the lowest frequency above noise. Assuming that the sample 
time is sufficient that the Fourier transform has converged, we define the threshold as 10% 
of the highest peak present in the power spectrum. For the simple graphs studied here 
the power spectra computed contain only a few frequencies and this functional definition is 
adequate. We expect, however, that for graphs with less symmetry the threshold may need 
to be modified. 

Using u s it is possible to determine the "period" of P s (t), that is the approximate integer 
time difference between minima in P s (t). To be of any use in the search context the search 
probability must reach a high value during the first period. We therefore quantify the success 
of the search as P max := max{P s (t) : 1 < t < 27i/u s ,t G Z}. 

The quantum search procedure resembles the wavelike propagation of a (phase inverted) 
perturbation over a graph. The perturbation originates from the marked vertex at each time 
step and results in time-dependent amplitude amplification at the marked vertex, measured 
by P s (t). The maxima in P s (t) occur when the probability amplitudes constructively inter- 
fere at the marked vertex, which is highly dependent on the structure of the graph. The 



amplitude and time dependence of P s (t) are therefore determined by both the local and 
global structure of the graph. In turn, these quantities may be used to provide information 
about this underlying structure. Indeed, using a similar procedure, Douglas & Wang [3] 
gave evidence that the information contained in these amplitudes was sufficient to distin- 
guish pairs of non-isomorphic graphs for all cases tested. It is reasonable therefore to expect 
that the centrality of a vertex should affect both P max and u s . 

III. CENTRALITY 

We now define the measure of centrality considered in this paper. The random walk 
centrality (RWC), introduced by Noh & Reiger (20}] is a measure of centrality designed to 
represent the relative speed with which a given vertex can receive and send information over 
a network. Pij(t) is defined as the probability of a classical random walker starting at the 
vertex Vi to be at the vertex Vj after a time t. The random walk centrality of a vertex Vj is 
then defined as 



where Tj = ^2tLo{Pjj{t) — Pj°} an d Pj° '■ = lim^oo P%j{t) which is the same for all i. Pij{t) 
is calculated using the master equation, 



Here A is the adjacency matrix of the graph and dk is the degree of vertex vj.. This definition 
was proposed in the context of complex networks and there is an implicit assumption that the 
expression for Tj converges, i.e. lim^oo Pij(t) exists. The master equation (Eq. [6]) describes 
a random walker that changes its location at each time step. Therefore, for graphs which 
contain closed walks of only even length, only sites at even distances from the start site 
are occupied after even times. In this case the probability distribution does not converge, 
which prevents the calculation of Tj (Eq. This difficulty is easily overcome by redefining 
the master equation in terms of the lazy random walk, such that the lazy random walker 
only moves at each time step with probability of 1/2. This leads to the following master 
equation, 



RWC, : 



poo 
3 



(5) 




(6) 
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FIG. 1. The third generation 3-Cayley tree. Vertices are ranked and labeled according to their 
random walk centrality. Non-labelled vertices are structurally equivalent to one of the labelled 
vertices. 



n this case, lim^oo Pij(t) = dj/N on all connected graphs regardless of the initial state 



21] . On the Sierpihski gasket, where both the lazy RWC and the normal RWC are defined, 
we find that they show the same qualitative behaviour. For the remainder of the paper, 
RWC will refer to the lazy random walk centrality. RWC is an example of a closeness 
centrality measure, in that larger values of RWC are associated with closeness to the centre 



of a network 



22] 



IV. TOPOLOGY AND CONNECTIVITY OF GRAPHS CONSIDERED 

We now describe the structures considered in this study. The n-th generation <i-Cayley 
tree (as shown in Fig. [1]) is a tree of n levels in which all vertices on the interior have 
degree d{ = d. The outermost layer of the tree is called the surface. All vertices on the 
surface are called leaves and have d% = 1. The total number of vertices in a Cayley tree is 
N = (d x (d — l) n — 2)/(d — 2). The structure is defined more rigorously in the Appendix. 
Quantum walks on Cayley trees, or the closely-related "glued-trees" graph of Childs et 



16 



23H26I]. and discrete-time 



22H2A 



al. have been extensively studied in continuous-time 
but not yet for discrete-time quantum search using marking operators. 

A first generation regular hyperbranched fractal (RHF) (as shown in Fig. EJ) of function- 
ality / is a star graph consisting of a central vertex connected through / edges to / surface 
vertices. To construct a second generation RHF, / copies of the first generation RHF are 
connected to the core first generation RHF through a single leaf-leaf edge. This procedure 



FIG. 2. (Color online) The third generation regular hyperbranched fractal of functionality, / = 3 
(RHFa^). The cul-de-sac vertices are highlighted. 




FIG. 3. The third generation joined 3-Cayley tree obtained from the equivalent Cayley tree by 
connecting surface vertices. Vertices are labelled as for the equivalent Cayley tree. Non-labelled 
vertices are structurally equivalent to one of the labelled vertices. 



is repeated n times for an n-th generation RHF. The number of vertices in an RHF there- 
fore grows exponentially with the generation and is given by the formula, N = (f + l) n . 
The maximum degree of a vertex in an RHF is /. Star graphs (RHFw) and other RHF n j 
have been studied previously for continuous-time quantum walks in [30|, l3l| . but not for 
discrete-time quantum walks. 

While the first generation Cayley trees and RHF! j are identical, the structures are dis- 
tinct for n > 1. An important difference in the context of this study is that for n > 2, RHFs 
contain so-called cul-de-sac vertices, which are leaves that are connected to main paths. 
Note that we do not consider the leaves of the Cayley tree to be cul-de-sac vertices since 
their neighbours do not lie on main paths. 
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FIG. 4. The N = 45 Husimi cactus obtained as a dual structure to the fourth generation 3-Cayley 
tree. Vertices are ranked and labeled according to their random walk centrality. Non-labelled 
vertices are structurally equivalent to one of the labelled vertices. 




FIG. 5. (Color online) The second generation Sierpihski gasket. Vertices are ranked and labeled 
according to their random walk centrality. Non-labelled vertices are structurally equivalent to one 
of the labelled vertices. Highlighted vertices are peripheral (least central as measured by RWC). 



We also study structures that contain more than one simple path between all pairs of ver- 
tices. Specifically, we consider the joined Cayley tree, the Husimi cactus and the Sierpihski 
gasket. The joined Cayley tree is obtained from the Cayley tree by adding edges between 
surface vertices (as shown in Fig. |3]). For trees with d = 3 this results in a 3-regular graph. 
The Husimi cactus (Fig. 0J is a dual structure to the Cayley tree, constructed by placing a 
vertex at each edge of the corresponding Cayley tree and connecting vertices that represent 
adjacent edges in the Cayley tree 3^]. The second generation Sierpihski gasket is shown in 
Fig. [5] and the structural details are explained in detail in 33j] . The number of vertices of 
an n-th generation Sierpihski gasket is N = |(3 n + 1). 
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FIG. 6. (Color online) a) Numerical results for the success probability P s (t) on the third generation 
3-Cayley tree (stucture shown in Fig. [1]) for a central marked vertex (vertex 1). b) P s (t) for a 
peripheral marked vertex (vertex 11) on the same graph. 

V. QUANTUM SEARCH ON CAYLEY TREES 



Firstly, we provide an example of the success probabilities obtained for two inequivalent 
vertices on a Cayley tree. We then present analytical results for the success probability 
P s {t) for the (i-Cayley tree when the central vertex is marked. Numerical results are then 
presented for P max and u s on various Cayley trees when a non-central vertex is marked. 
These results are compared with the centrality of the marked vertex. For all simulations in 
sections [V] ED and IVIIl the initial state is given by Eq. |2j 



A. Quantum walk-based search characteristics 

We begin with an example of the success probability obtained on the third generation 



3-Cayley tree (structure shown in Fig. [p. Fig. 6(a) shows the success probability P s (t) as a 
function of time when the central vertex is marked. Fig. 6(b) shows P s (t) for a peripheral 
marked vertex on the same graph. Comparing these plots, it can be seen that P 8 (t) is 
quasi-periodic in both cases, with a smaller "period" when the central vertex is marked 
in comparison with a peripheral marked vertex. P s (t) also has a greater maximum and 
average amplitude for the central marked vertex. In the following, we examine how the 
lowest frequency u s and the maximum amplitude P max vary on Cayley trees. 
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FIG. 7. (Color online) a) Analytically obtained curves for the search frequency ui s (n,d) for n = 2 
and 3 on a Cayley tree with a central marked vertex. Also plotted are numerical data obtained 
by direct simulation of discrete-time quantum search, b) Random walk centrality for the central 
vertex on the same Cayley trees. 

B. Central marked vertex 

In the Appendix we show that a quantum walk on an n-th generation (i-Cayley tree can 
be mapped to a one-dimensional quantum walk. We then study a finite one- dimensional 
walk with two reflecting boundaries to derive the following expression for u s , 



2Vd n - 1 - 1 
1) 



u s = u s (n, d) := arctan ( ) (8) 



for n — 2, 3 and 3 < d < oo. 



This expression is plotted in Fig. 7(a) together with our numerical results obtained by 
direct simulation of discrete-time quantum search. The numerical and analytical results are 
in perfect agreement. Finding analytical solutions for uj s becomes more difficult for n > 4 
as large matrices must be diagonalised. For the second generation Cayley tree (n = 2), we 
are able to derive (see Appendix) the following expression for P s (t), valid for 3 < d < oo, 

1(1 W))|2 = 4(i + rf2 ) ( 1 + d 2 + (d- l) 2 cos(7rf) 

- (d 2 - 1) cos(w s t) + 2^Jd(d- l)sin(w s t) 

- (d 2 - 2d - 1) cos[(tt - u a )t] 



- 2y / d{d - 1) sin[(7r - u) a )t]}. (9) 
As can be seen from Eq. [91 P s (t) contains only three frequencies, 7T, u s and n — oj s . It is 
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FIG. 8. (Color online) a) Numerical results for the search frequency co s on a Cayley tree with 
central marked vertex. oj s is plotted as a function of the number of generations n, for d = 3 
(solid/circles) and d = 4 (dashed/squares), b) Random walk centrality for the central vertex on 
the same Cayley trees. 



interesting to note that u; s (2, 4) = 7r/3, which means the only frequencies present in P s (t) 
are {vr/3, 27r/3, 7r}. Thus P s (t) is exactly periodic with period 6. 



The search frequencies oo s (2, d) and w s (3, rf) are plotted in Fig. 7(a) Although for fixed n 
the distance between the central node and the surface is fixed, it can be seen from Fig. 7(a) 
that U) s (n,d) decreases monotonically with increasing branching rate. According to EaJSJ 
for large d, u s tends toward 0. This is consistent with the results of Carneiro et al. [29], 



who generalised a discrete-time formulation of Childs' glued-trees algorithm [23] to include 
arbitrary branching rate d and found that as d — > oo, a quantum walker initially localised at 
the central node of the tree oscillates between the central node and the first level of the tree. 
While it is not strictly valid to compare absolute values of RWC between different graphs, 
we find numerically that RWC of the central vertex of a Cayley tree is also a monotonically 



decreasing function of d for fixed n, as shown in Fig. 7(b) In this sense, u s shows the same 



behaviour as RWC for the central marked vertex of a Cayley tree. 



For n > 4 we only have numerical results for u s , which are shown in Fig. 8(a) The 
search frequency u s on a Cayley tree with a central marked vertex decreases as the number 
of generations in the tree increases. Thinking of the search procedure as generating a phase- 
inverted perturbation at the centre of the graph at each time step, these perturbations must 
constructively interfere at the centre of the graph in order to produce a maximum in P s {t). 
To constructively interfere, these perturbations must be reflected at the surface and thus 
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FIG. 9. (Color online) Analytical results for P ma , x on a second generation Cayley tree with a 
central marked node as a function of branching rate d (solid line). Also plotted are numerical 
results (circles) obtained by direct simulation of quantum search for 3 < d < 10. 



the time between maxima in P s (t) depends on the distance between the central vertex and 
the surface of the graph, which increases with the number of generations, n. As shown in 



Fig. 8(b) , RWC of the central vertex decreases monotonically with n for fixed d and thus 
displays the same behaviour as u s for the Cayley tree. 

We now study how P max varies with d for a central marked vertex on the second generation 
Cayley tree. By evaluating Eq. [9] at integers t ~ n/u s (2,d), we can find -P max for arbitrary d. 
The results are shown in Fig. [9] for 3 < d < 40. Numerical results agree with our analytical 
results for 3 < d < 10. Fig. [9] shows that P max has a much more complex dependence on 
d than was seen for u s . This complex dependence on d arises because the frequencies it, 
u s (d) and tt — u s (d) present in P s (t) do not necessarily constructively interfere within the 
first period for all d. While RWC and u s both decrease with increasing branching rate, we 
see that -P max generally becomes larger as d increases. This can be understood as arising 
from the greater localisation of the quantum walk close the central vertex as branching rate 
increases. 



C. Non-central marked vertices 

The analysis in the Appendix is dependent on being able to map the quantum walk on 
a Cayley tree to a walk on a line. When a non-central vertex is marked, the same mapping 
cannot be used and we do not have analytical results. 
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FIG. 10. (Color online) Numerical results for P max , uj s and RWC on the third generation 3-Cayley 
tree (all data normalised for comparison). Data plotted as a function of vertex number, V{. 



Instead we present numerical results for P max and oo s on Cayley trees where the marked 
vertex is not necessarily central. P s {t) was computed via direct application of U to the 
probability amplitudes of the subnode states \vi, Cj) on which the quantum walk takes place. 
Numerical methods were then used to compute the discrete Fourier transform and u s was 
obtained as defined in Sec [TTJ For a given graph, this was repeated for all possible positions 
of the marked vertex. Fig. [10] shows the dependence of P max and u s on the centrality of 
the marked vertex for the third generation Cayley tree with d = 3 (structure shown in 
Fig. CE]). As can be seen in Fig. (TOj a more central vertex on Cayley tree gives rise to a 
success probability P s {t), which has a larger minimum frequency, u s , and attains a greater 
maximum amplitude, P max - The results for all other Cayley trees studied were analogous. 

On Cayley trees, RWC becomes smaller for vertices which are further from the central 
vertex because a random walker starting from a less central vertex takes, on average, longer 
to visit all vertices. We studied numerically the following Cayley trees: {d = 3,n G [2,8]}, 
{d G [4, 5], n G [2, 5]} and {d G [6, 10], n G [2, 3]}. For all Cayley trees tested we found that 
-Pmax, (P s ) and uj s decrease monotonically with RWC. 

These results are consistent with those observed by Mulken et al. j^] who modelled 
exciton transport on the Cayley tree by continuous-time quantum walk and found that 
excitations which were initially centrally located propagated throughout the structure much 
more rapidly than excitations which were initially peripheral. 
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FIG. 11. A comparison of the structures of a) the second generation 3-Cayley tree and b) RHF2,3. 
Vertices are ranked and labelled in order of decreasing random walk centrality. 
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FIG. 12. (Color online) Numerical results for -P max , w s , and (P s ) on RHF2,3. All data is normalised 
for comparison and vertices V{ are ordered with decreasing random walk centrality. 

VI. QUANTUM SEARCH ON REGULAR HYPERBRANCHED FRACTALS 

We now study RHF in order to investigate how greater structural complexity affect both 
P max and uj s . 

Firstly, we consider the simple case of n = 2, / = 3 (RHF^s). As shown in Fig. [HI the 
only difference between this fractal and the 2nd generation 3-Cayley tree is extra vertices 
of degree 2 on each branch, between the central vertex and level 1. In this case, like on 
Cayley trees, we find that u s decreases monotonically with RWC. However, as shown in 
Fig- EE21 neither P max nor (P s ) decrease with centrality. In fact there are non-central vertices 
in RHF 2 3 at which P max and (P s ) attain a greater value than at the central vertex. This 
implies that P max and (P s ) are not determined by the centrality of the vertex. It is remarkable 
to consider that the addition of these 6 vertices symmetrically about the central vertex can 
perturb the quantum walk dynamics to such an extent that the the relationship between 
centrality and P max seen on Cayley trees is no longer maintained. 
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FIG. 13. (Color online) Numerical results for P majX , u s and RWC on RHF 33 . All data is normalised 
for comparison and vertices Vi are ordered with decreasing random walk centrality. The (*) denotes 
sets of cul-de-sac vertices. 

Fig. [T3l shows the results of the calculations forn = 3, / = 3 (RHF 3j3 - structure shown in 
Fig. [2]). These results show that the overall trend is again decreasing uj s with RWC but that 
the decrease is not monotonic. This implies that the search frequency on RHF is not solely 
determined by the centrality for some vertices. The vertices which break the trend have a 
particular structure, they are vertices which are connected by a single edge to a vertex which 
lies on a main path of the graph, so-called cul-de-sac vertices (highlighted in Fig. [2]). The 
much smaller than expected values for u s at these cul-de-sac vertices is the first evidence 
that the local structure surrounding the marked vertex can greatly affect the frequency of 
the success probability. We now investigate this further. 

We construct the graph RHFjj'g by adding a single extra vertex to RHF 33 attached 
through a single edge to vertex 16, which is one of the cul-de-sac vertices in RHF 3;3 (see 
Fig. E]). The calculated results of P max , co s and RWC for RHF 3j3 and RHF33 are shown in 
Table [H Note that in RHF 33 , vertices 14, 15 and 16 are structurally equivalent and have 
the same values for P max , u s and RWC. The additional vertex attached to vertex 16 in 
RHF^ 3 results in decreasing RWC for that vertex. Given that RWC is calculated using the 

16 





RHF 33 


RHF+3 


change 


Vi = 14, 15 


p 

1 max 

U s 

RWC 


2.065 x KT 1 
4.533 x 10~ 2 
1.166 x 1(T 3 


2.023 x 10" 1 
4.574 x 10~ 2 
1.147 x 10~ 3 


-2.0% 
+0.9% 
-1.7% 


Vi = 16 


p 

1 max 

RWC 


2.065 x 10" 1 
4.533 x 10- 2 
1.166 x 10~ 3 


2.077 x 10- 1 
4.533 x 10~ 2 
1.106 x 10~ 3 


+0.6% 
0.0% 
-5.4% 



TABLE I. Changes in the values of -P max > w s an d RWC between RHF3 3 and RHF^g, which has 
an extra vertex added to vertex 16 (see Fig. [2] for structure of RHFs^). 

relaxation time, this is expected. Since centrality is a relative measure among vertices, it is 
more informative to compare values within the same graph. We therefore compare values for 
vertex 16 with those for vertices 14 and 15, which are all equivalent in RHF3 3. For RHF^g, 
we compute the ratios, 



p(16) 

-'max _ 1 n9 
(14,15) ~ X - UZ ' 



-Pmax 

t ; (16) 
, ,(14,15) 

RWC (16) 
RWC (14 - 15) 



0.99, 
0.96, 



and we see that in real terms, the effect of the node addition to vertex 16 is to decrease both 
u s and RWC. So we see again that u s appears to behave like RWC. 

These data suggest that, on symmetric trees, u s is largely determined by RWC but the 
local structure surrounding the marked vertex can have a dramatic effect. It should be 
noted that all vertices on RHF which are observed to break the trend of decreasing u s with 
centrality have degree 1. 
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FIG. 14. (Color online) Numerical results for P ma , x , oo s and RWC on the Husimi cactus on N = 45 
vertices. All data is normalised for comparison and vertices Vi are ordered with decreasing RWC. 

VII. QUANTUM SEARCH ON STRUCTURES WITH SIMPLE CYCLES 

We now move away from trees to consider quantum walk-based search on structures which 
contain simple cycles. In particular, we consider the Husimi cactus, the joined Cayley tree 
and a Sierpihski gasket. 

While a tree contains no simple cycles, a cactus is a graph in which any edge belongs 
to at most one simple cycle. In this sense, cacti are the most "tree-like" graphs which 
contain simple cycles. On Husimi cacti derived from 3-Cayley trees, all simple cycles have 
length = 3. Studies of continuous-time quantum walks on these Husimi cacti have found 



that they display similar dynamics to Cayley trees [321 . Fig. [TH shows u s , P max and RWC 
for the N = 45 Husimi cactus (structure shown in Fig. H]). Like on the Cayley tree, it can 
be seen that the search frequency u s decreases with the centrality of the marked vertex. 
However, in contrast to the Cayley tree, P max does not monotonically decrease with RWC. 
The results for the N = 21 Husimi cactus are analogous. 

Fig. [151 shows u s , P max and RWC for the joined Cayley tree of generation 3, with d = 3. 
We again see that u s and P max generally decrease with centrality in a similar manner to 
the unjoined Cayley tree. In the joined graph (Fig. [3]), vertices 11 and 12 are no longer 
structurally equivalent since vertex 11 lies in simple cycles of length = {3, 7} and vertex 12 
lies in simple cycles of length = {3, 5}. It can be seen from Fig. [15]that the search procedure 
for marked vertices 11 and 12 now produces different values for RWC, P max and u a . It should 
be noted that, for these vertices, P max and RWC display the opposite trend to cu s . This is 
different to the correlation observed on trees and Husimi cacti, where we saw that uj s was 
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Vi 



FIG. 15. (Color online) Numerical results for P max , uj s and RWC on the third generation joined 
3-Cayley tree (all data normalised for comparison). Data plotted as a function of vertex number, 



Vi, which is derived from the numbering of the unjoined Cayley tree. Note that vertices 11 and 12 



are no longer equivalent in the joined Cayley tree. 




0.0 

2 4 6 8 10 12 14 



FIG. 16. (Color online) Numerical results for P max , uj s and RWC on the second generation 
Sierpihski gasket. All data is normalised for comparison and vertices V{ are ordered with decreasing 
RWC. 



usually highly correlated with RWC. 

The results in Fig. [16] for the Sierpihski gasket show a much more complex relationship 
between P max , u s and centrality to that seen on the other structures studied. We expect 
that the presence of simple cycles of various lengths in the structure results in interesting 
interference effects for quantum walks on these graphs. This interference produces much 
more complex behaviour of P s (t), which makes any relationship between P s (t) and the 
centrality of the marked vertex harder to observe. It should be noted from Fig. [16] that the 
peripheral (least central) vertices do not produce the smallest values of P max or oj s in this 
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case. 



VIII. DISCUSSION & CONCLUSIONS 

We have demonstrated that the quantum search procedure for a marked vertex on a 
Cayley tree results in a time- dependent success probability P s (t) that attains a greater 
maximum value and has a higher minimum frequency for marked vertices which are more 
central. Our study was extended to regular hyperbranched fractals where it was found that 
the maximum value of P s {t) does not necessarily decrease with centrality. We also found 
that the minimum frequency of the search probability u s on RHF is strongly related to 
centrality and that exceptions to this trend are caused by the local structure surrounding 
the marked vertex. We therefore conclude that the success probability P s (t) for a marked 
vertex on a highly symmetric tree contains information about the global structure of the tree 
and the overall position of the marked vertex within the tree. It also contains information 
regarding the local structure surrounding the marked vertex. 

We say that two vertices are structurally equivalent if they have the same structural 
relationship to all of the other vertices of the graph (e.g. on a Cayley tree, all vertices 
which are a given distance from the centre are equivalent). We then note that all equivalent 
vertices on the graphs studied in this work have the same values of P max and u s and all 
inequivalent vertices have different values, suggesting that P max and u s could be used to 
partition the vertex set into structural equivalence classes. It should be noted however that 
this partitioning would not be possible for all graphs, since strongly regular graphs contain 
inequivalent vertices that produce identical success probabilities. These graphs therefore 
give identical values for u s and P max for inequivalent vertices. We find that for all vertices 
of degree > 1 on the trees and cacti studied here, u s can be used to order the vertices with 
decreasing centrality, as measured by RWC. 

The random walk centrality considered in this paper is an example of a closeness cen- 
trality measure, i.e. it measures the closeness of a vertex to the centre of the graph. There 

□ 

are, however, other classes of centrality measures [22]. One of these is betweeness centrality, 
which represents the capacity of a vertex in a network to withhold information if it were 



removed from the network 221 ] . While betweeness and closeness are equivalent on a Cayley 



tree, this is not true for RHF. As seen in the studies on RHF, the values of u s obtained for 
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the cul-de-sac vertices were unusually low when compared with the corresponding values of 
RWC. This could be related to the low betweeness centrality of these vertices. It would be an 
interesting subject for further study to investigate the relationship between characteristics 
of P s (t) and other measures of centrality 

It is also worth noting that P s (t) can be classically computed with time complexity 
0(tN 2 log N) where t is the simulation time and N is the number of vertices. Therefore, u s 
can be classically computed with time complexity O (t 2 log t N 2 log N) . Numerically, we find 
that for a 1% error in u s for typical values of u s on a 3-Cayley tree, t « 0(N°' 7 ) and thus 
the time complexity becomes 0(iV (log N) 2 ). This is comparable with the complexity of 
the equivalent RWC calculation which scales as 0(tN 2 ) 0(N 3A ). 

It is evident that a tremendous amount of information about the structure of the graph, 
and the overall position of the marked vertex is present in P s (t) and we suggest that the 
"lowest frequency" considered here is perhaps the simplest example of extracting this in- 
formation. We propose that it would be extremely interesting to consider more detailed 
Fourier analysis of P s (t) to uncover the origin of the other frequencies present and hence 
determine what other information is obtainable. It would also be interesting to determine if 
the centrality information, which was easy to obtain via u s in the case of simple symmetric 
trees, can be extracted for more complex cases. 
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Appendix: Analytical results for oj s on Cayley trees. 

We map the quantum walk on a CT to a finite one- dimensional quantum walk and derive 
an expression for u s when the central vertex is marked for n — 2, 3. The first part of the 



analysis follows the method of Chisaki et al. 27] however the time-evolution operator used 



here is modified to incorporate the marking operator and reflection at the surface vertices 
(Chisaki et al. consider an infinite tree). The initial state is also different, therefore the 
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|£i£o> |£2£o) |£q£i) |£2£i> |eo£a> |£i£2> 



FIG. 17. The second generation 3-Cayley Tree, shown with the vertex states \g) G Tip and the 
subnode states \g,e) G %-p <8> He- The shaded regions enclose groups of subnode states that make 
up { I a;) ou t, |^}m}; shown on the right. 

complete derivation is presented. 

1. Defining the discrete-time quantum walk on a Cayley tree 

An n-th generation Cayley tree of order d is an undirected graph in which every vertex 
is connected to d others except for vertices at a distance of n from the centre. The vertex 
set, V is defined as the group of elements generated by the set £ = {e , . . . , Sd-i\ with the 
constraint £j 2 = 1, i — 0, . . . , d — 1, 

V = {e^Ei^ . . . e i2 e h : < k < n, G E, 
and ij + i ^ ij for j — 1, 2, . . . , k — 1}. 

The reduced word length, \g\ of a vertex g is the number of elements used when g is written 
as a reduced product of £j G £ (e.g. I^^i^o^l = Iki^o^l = l^i^o^l = 3). Using this 
construction, the vertices at level k in the tree have \g\ = k. Vertices g and h are connected 
if and only if gh~ l G £. This is illustrated in Fig. [TTJ 

The quantum walk takes place on the subnodes of the graph which belong to the Hilbert 
space Tip <S> He- We write these states {\g, Sj) : g G V,Ej G £}. In this basis the action of 
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the unitary Grover coin and shift operators which drive the quantum walk can be written 

as, 



[l®C)\g,e) 



ErGs(~ tf er + 2/d)|p,r) , \g\ < n 



S\g,s) = \eg,e). 



\9\ = n 



(A.l) 
(A.2) 



The coin operator (Eq. IA.1I) applies a Grover coin to all vertices except those at the surface 
which are left unchanged, while the shift operator (Eq. IA.2j) swaps the probability amplitudes 
between connected subnodes. With the perturbed time evolution operator U' := 5(1 ® C), 
one step of the discrete-time quantum walk on the Cayley tree becomes, 



U'\g,e) = { 



J2r^(-^r + 2/d)\rg,r) ,\g\<n,g£M 
\eg,e) , \g\ = n,g$ M 

-\eg,e) , g G M. 



(A.3) 



2. Mapping the quantum walk on a Cayley tree to a quantum walk on a finite line 



As described in 23j, [24|, |27|, l22| , in certain cases, the symmetry of the initial state allows 
us to map the quantum walk on a Cayley tree to a walk on a line. Define the sets 

E + {x) = {(g, e) e V x E : \g\ = x, \eg\ = x + 1}, 
E ~( x ) = {(g,e) e V x S : \g\ = x, \eg\ = x - 1}, 

and consider the subspace W C % spanned by following states, 

1 



|>^)out 
|-^)in 



< x < n — 1 



1 < x < n. 



(g,e)eE-(x) 

as shown in Fig. [TT1 If the initial state |^o) °f t ne quantum walk on a Cayley tree is an 
equal superposition of all vertex states, equally divided among the subnodes of each vertex, 



l*o> 



d-2 



(dx(d- 1)" -2) 



E E 

g&V e:\eg\<n 



Vdjg) 



(A.4) 
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where d(g) is the degree of the vertex, then |\I/o) can be written as a superposition of 

\ 1 •£) out; |-^)in}- 

n-1 



x=0 

n-1 



x=l ' 



We now specialise to the case where there is a single marked vertex at the centre of the 
tree (M = {1}) and consider the action of U' on the subspace %' C H-p <S> He- 
For 1 < x < n, 



U'\x) in = (l-l)\ x - l) out + 2v ^^ |x + l) in , (A.6) 

U'\x) ont =-(L-l)\ x+ l) in + 2y/T ~ l \x - l) out , (A.7) 
a a 

and on the boundaries, 

U'\0) out = -|l) in , (A.8) 
U'\n)- m = \n-l) ont . (A.9) 

Note that |* ) eW ^ = (E/')*l*o) g H', Vt g Z+. 

From this subspace the quantum walk on the tree can be mapped to a quantum walk on 
a finite line with reflecting boundaries. Define the Hilbert space H — {\x, A) : x G Z + , A G 
{L,i?}} and the subspace W C H spanned by {|0, L), |1, R), |1, L), . . . , \n — 1, R), \n — 
l,L), \n,R)}. Now consider the 1-1 association %' -H- defined by, 

\x) ont -H- \x,L), \x)- m -H- (A. 10) 

We also map the operator [/' <H- [/' so that the action of U' on the states {|x) in , |x) out } is 
equivalent to the action of U' on the states {\x, R), \x, L)}. The quantum walk on a Cayley 
tree subject to the above conditions is then seen to be equivalent to a finite one-dimensional 
quantum walk with a biased coin and asymmetric initial distribution. 

Using the mapping IA.10| we can express IA.6[ IA.7[ IA.8[ IA.9I in the new basis by defining 
W := S(t <g> C), where C = C(x) is defined by, 

(1 ®C)\x,A) = \x) ®H(x)\A). (A.ll) 
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where H(x) 



, < x < n , 
, x = n 



2Jd^T 



h 



- 1 



d d 
'2 i\ 2yfd^l 



■(1-1) 



and <j\ 




S\x,A) 



(A.12) 



S acts on W as follows, 

+ A = i? 

\x-l,L), A = L 

It can thus be seen that the quantum walk on the Cayley tree with a central marked 
vertex starting with equal probability at all vertices is equivalent to a quantum walk on a 
finite line with a perfectly reflecting boundary at x = n and a ix phase shift upon reflection 
at x = 0. The reflections are acheived by careful choice of the coin operator. 



3. Analytical results for P s (t) and ui s for the 2nd generation Cayley Tree 

Consider the 2nd generation <i-Cayley Tree with a central marked vertex. Let B be a 
matrix representation of U' on the complete, ortho normal basis {|0, L), \1,R), \1,L), |2, R)} 
of H'. Let G {|0, L), |1, R), |1, L), |2, i?)}. 



then (Bij) = [(<t>i\U'\<l>j) 



(° 


2 1 2 V / 5^T 
d d 


o\ 


-1 














1 


u 


2^/d^T I 2 
d 1 d 


V 



where the states are ordered as above. 

We would like to find a closed form expression for \^(t)} = ([/')* |^ )- In particular we are 
interested in the probability of finding the walker at the central marked vertex as a function 
of time. 
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which can be written in matrix form as, 



d(t) = (BitfcjiO), (A.13) 

where Cj(i) := (fal^f (t)) are column vectors. Since U' is a unitary operator and the basis 
\4>i) is orthonormal, it follows that B is a 4 x 4 unitary matrix and therefore has 4 distinct 
eigenvalues of unit norm corresponding to 4 orthonormal eigenvectors js4]- Furthermore, 
B is unitary similar to the diagonal matrix D, which contains the eigenvalues of B as its 
diagonal elements. Let P be the unitary matrix containing the eigenvectors of B as its 
columns. From Eq. flA.13j) and the unitarity of P, 



P\ iCl (t) = P^kiiBijYPjiPhjCjiO), 

phMt) = (D kl yp\ jCj (o). 

Now writing P\iCi{t) =: Vk(t) we have, 

v k {t) = p w )V0). 



(A.14) 



But D = diag(Afc) where are the eigenvalues of B. So Dm = for k ^ I =>- Vk(t) — 
(DkkYv k (0) or equivalently, 

v k (t) = X k t v k {0). (A.15) 

The eigenvalues of B are, 



X k = ±e ±n/2 , where Q = arct 



an 



(A.16) 



d-2 

We can now solve for the time evolution of Vk(t) for k — 1, . . . , 4 and use P to find c&(t) := 
{4>k\ty(t)), from which we obtain The eigenvalues are found by diagonalising B, 



i.e. 



D = P^BP = dia K (-e- lQ / 2 , -e lQ / 2 , e^ 2 , e^ 2 



(A.17) 



where the unitary change of basis matrix P contains the orthonormal eigenvectors of B as 
its columns. 

I i —i —i % \ 



p= 1 - 

2 



_ e iQ/2 _ e -in/2 g -»n/2 gift/ 2 

V 1 1 1 1 / 



(A.18) 
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The probability amplitude at the central (marked) vertex is thus 

(l|tt(t)> EE (0,L|*(t)> = Cl (t) =J2PlMt) 

i 

-*n/2)t Vl ( ) _ i(_ e * n /2)* T , 2 (o) 



i — e 



ze^^O) + ie-*"*/ 2 ^) 



where fj(0) represent the initial state in the eigenvector basis and are found using Eq. IA.5I 
and the unitary transformation P'. Taking the absolute value squared, 

|(l|^(t))| 2 =^{|v| 2 + 2 Re(«iv 4 +v 2 v 3 ) cos(Trt) 

— 2 Re(f 1W2 + ^3^4) cos(fit) 

— 2 Im(-uif2 + sin(fii) 

— 2 Re(f 1^3 + ^2^4) cos[(7r — fl)t] 
+ 2 Im^i^ + -^2^4) sin[(7r — 

We now see that the lowest frequency in the success probability, u s = Q. The initial condition 
(Eq. IA.5I) can be converted to the eigenvector basis and the expression (Eq. IA.16I) for Q(d) 
can then be used to solve for t>j. Upon simplification, we obtain Eq. [9] (Sec. IV Bp . valid for 

3 < d < 00. 



4. Analytical results for P s {t) and ui s for the 3rd generation Cayley Tree 

We perform the same analysis for the case n = 3 where this time U' induces a 6 x 6 
matrix B. We find that the eigenvalues of B are, 

A fc = ±i, ±e ±n/2 , where Q = arctan f ^~^~ ) ■ ( A ' 19 ) 

As in the previous case, the eigenvalues of B uniquely determine the frequencies present in 
P s (t). Since i = e i7r/2 and Q < ir/2 for 3 < d < 00. It follows that u s = Q. 
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